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Abstract 
We  study  the  question  of  when  the  power  spectrum  of  a  stationary- 
random  process  is  uniquely  determined  by  the  values  of  the  process 
measured  in  various  ways  at  discrete  instants  of  time.     As  is  well 
known,   in  the  case  of  equi-spaced  sampling,   the  power  spectrum  is 
not  \iniquely  determined  by  the  sampled  data,   i.e.  aliasing  is  present. 
After  a  review  of  the  equi-spaced  case,  we  examine  the  case  where 
the  time  markers   are  not  equi-spaced,  but  are  subject  to  jitter,  and 
find  that  in  general  aliasing  persists.     We  then  consider  the  case 
where  the  time  markers  are  generated  by  an  additive  scheme,  whereby 
each  sampling  time  is  derived  from  the  previous  one  by  the  addition  of 
an  independent  random  variable,  with  characteristic  function  9^(3).     It 
is  found  that  additive  random  sampling  is  alias-free  if  0{s)  is 
one-to-one  on  the  real  axis,  but  not  alias-free  if  ^{s)  takes  the 
same  value  at  two  distinct  points  of  the  open  upper  half -plane. 
Various  alias-free  sampling  methods  are  exhibited,    notably  Poisson 
sampling. 
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1,  Introduction, 

A  problem  of  great  practical  interest  is  the  determination  of 
the  power  spectrum  fCco)  of  a  stationary  random  process  x(t)  from 
measurements  of  x(t)  made  at  discrete  instants  of  time.  As  is  well 
known  (see  e.g.  [l] )  the  usual  method  of  sampling  at  equally  spaced 
instants  of  time  t  =  nh,  n  =  . . .,  -2,  -1,  0,  1,  2,  .. . ,  does  not 
permit  unambiguous  detemdnation  of  F(co),  unless  F(a))  is  known  in 
advance  to  lie  in  the  Nyquist  band  N,  defined  by  |co|  <  rt/h;  otherwise, 
a  whole  class  of  different  pov;er  spectra  are  compatible  with  the  sampled 
values  of  x(t).  In  the  usual  terminology,  these  different  power 
spectra  are  called  aliases  of  F(co)  .  Thus,  it  is  clear  that  if  we 
hope  to  find  alias -free  sampling  methods,  i.e.  methods  which  lead  to 
unambiguous  determination  of  F(cij),  we  must  sample  at  unequally 
spaced  instants  of  time.  The  case  where  these  unequally  spaced  times 
are  generated  purely  arithmetically  seems  to  lead  to  great  analytical 
difficulties  .  However,  as  shown  in  the  present  paper,  the  case  where 
the  sampling  times  are  chosen  randomly  is  quite  tractable.  In  fact, 
we  shall  see  that  some  random  sampling  schemes  succeed  in  eliminating 
aliasing,  while  others  do  not.  In  particular,  the  important  case  of 
Poisson  sam.pling,  where  the  sampling  times  are  taken  to  be  the 
occurrence  times  of  the  events  of  some  Poisson  process  (e.g.  the  arrival 
times  of  shots  in  shot  noise  or  the  particle  emission  times  in  a 
radioactive  decay  process)  is  found  to  be  alias-free. 


Equivalently,  we  call  the  corresponding  correlation  functions  aliases 
of  one  another. 

See  remark  on  p.  219  of  [IJ , 
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2r  Periodic  Sampling. 

In  what  follows  we  shall  be  concerned  exclusively  with 
second  order  (or  wide  sense)  properties  of  x(t),  which  may  therefore 
be  taken  to  be  a  stationary  Gaussian  process.  We  assume  that  x(t) 
is  real,  with  mean  zero  and  continuous  correlation  function  C('c),  i.e. 

Ex(t)  =0  ,  Ex(t)x(t')  =  C(t-t')    , 

where  E  denotes  the  expectation  or  ensemble  average   .     3y  the  usual 
argument,  C(-r)  =  C(r)  and   lG(-t')l   <  C(0).     We  write  the  Wiener-Khint chine 
relations  in  the  form 

00  /OO 

(1)  C(t:)   =         exp(ioir)F((o)do     ,     F(ai)   =  ■—         exp(-ixof)C(^)df , 

•^00  '^CO 

where  the  even,  non-negative  function  F(oj)  is  the  power  spectrum  (more 

accurately,   the  power  spectral  density)  of  x(t).     It  will  be  assiuned 

(^  2 

that  F(co)  is  not  only  integrable,   as  follows  from        F(6i)doii  =  Ex  (t)   , 

-00 

but  also  square-integrable.  Symbolically  ,  we  have  F(co)^L  DL  .  It 

follows  by  Plancherel's  theorem  that  C('C)£'L  ,  and  that  the  second 

/R 

integral  in  (l)  may  have  to  be  interpreted  as     p*"^^* 


We  shall  always  assume  that  x(t;  has  the  ergodicity  needed  to   replace 
the  operation  E  by  suitable  infinite  time  averages.     In  the  case  of 
Gaussian  x(t),   it  follows  by  a  theorem  of  Maruyama   [2,  Thm.   h]  that 
x(t)  is  ergodic,   since  the  required  continuity  of  the  spectral  distri- 
bution function  of  x(t)   is  implicit  in  the  representation  (1).     In 
this  first  study  of  alias-free  sampling,  we  shall  not  study  the  variance 
of  finite-time  power  spectrum  estimates   (as  done  e.g.  in   [l]). 

By  f(co)^L^(a,b),   p  =  1,2,  we  mean  as  usual  that  the   (Lebesgue)  integral 
J    |f(co)pd£0  is  finite.     We  abbreviate  lP(-oo,oo)  to  L^, 
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Suppose  now  that  we  sample  x(t)  periodically  (i.e.  with  equal 
spacing)  at  the  points  t  =  nh,  n  =  ...,  -2,  -1,  0,  1,  2,,,.,  where 
h  >  0  is  the  sampling  interval.  This  generates  the  stationary  random 
sequence  x(t  )  =  x(nh)  ,  with  mean  zero  and  correlation  function 
(equivalently,  correlation  sequence) 

c^(n)  =  C(nh)   ,    n  =  ...,  -2,  -1,  0,  1,  2,... 

Since  0(1:')  is  even,  we  have  c,  (-n)  =  c,  (n).  The  discrete  analogs  of 

h       h 

(l), which  will  be  needed  later,  are 

/  1   °° 

(2)    Cj^(n)  =  J  exp(ixon)fj^(co)dco  ,   fj^(co)  =  ^  y~     exp(-iMn)c^(n)  , 

-n  n=-oo 

where  f,  (co),  the  power  spectinm  of  x(t  ),  is  defined  only  in  the 
interval  (-TT,n), 

We  now  give  two  different  ways  of  constructing  aliases  for  tlie 
case  of  periodic  sampling.  The  first  method  relies  on  the  fact  that 
any  two  correlation  functions  which  agree  on  the  "  lattice"  of  points 
t  =  nh  are  obviously  aliases  with  respect  to  periodic  sampling  with 

spacing  h.  To  construct  such  a  pair  of  correlation  functions,  let  h(f) 

1    2 

denote  any  real  even  function  in  L  f]  L  ,   possessing  a  continuous 

second  derivative  and  vanishing  at  the  points  t  =  nh  .  Then  the 


We  exclude  the  trivial  case  h(f )  £  0, 
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Fourier  transfoiri 

,00 


H(co)   =  ^ 


exp(-icot)h('tr)d'G 


/ 


-co 

no  O 

is  real,   even  and  in  L  0  L     (since  H(co)  =  Odwl"*  )  for  large  w).     If 
now  we  write 


(3) 


and 


a) 


H^(o3)   =  H(co)      ,     H(co)  >  0        , 
H_^(co)  =  0  ,     H(co)  <  0        , 


H_(co)   =  -H(w)   ,     H(w)  <  0        , 
H_(co)   =  0  ,     H(co)  >  0        , 


then  H_j_(a))  and  H__(w)  are  a  pair  of  even,  non-negative  functions  in  L  /|L  , 
whose  Fourier  transforms  C^{t)  and  C_(T;)  agiree  on  the  lattice  t     =  nh, 
and  hence  are  aliases  of  each  other.  '     Note  that  H^(co)  and  H  (co)  are 

not  only  distinct,  but  even  non-overlapping.     However,   this  method  can 
also  be  used  to  construct  aliases  which  lie  in  the  same  frequency  band. 

We  merely  note  that  all  the  functions 

(5)  H  (<o)  =  aH.(w)  +   (l-a)H  (co)         ,     0     <     a     <     1     , 

are  aliases  of  one  another.  In  fact,  more  generally,  any  element  of 


By  construction,  C_j_(f)  and  C^Cf)  are  automatically  positive  definite, 
and  hence  correlation  functions. 
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the  convex  hull  of  a  set  of  alias  power  spectra  is  another  alxas. 

Using  the  same  construction,  we  can  show  that  knowing  the 
correlation  functions  of  the  finite  family  of  random  sequences 
x(nh_ ),., .,x(nh  )  is  not  enough  to  eliminate  aliases,  i.e.,  there 
exist  pairs  of  correlation  functions  which  agree  on  any  finite  number 
of  lattices  of  equally  spaced  points.  It  is  merely  necessary  to 

choose  the  function  h(T)  in  the  above  constriction  to  be  any  real 

1   ? 
even  function  in  L  DL  ,  possessing  a  continuous  second  derivative 

and  vanishing  on  all  the  lattices  nh,,...,nh  ,  and  then  construct 

H  (w)  and  H  (co)  as  before.  Again,  more  aliases,  with  overlapping 

spectra,  can  be  constructed  by  using  {$) , 

The  second  method  of  constructing  aliases  is  based  on  a  detailed 

examination  of  the  relation  between  the  power  spectrum  fy,(")  of  "t^© 

random  sequence  x(nh)  and  the  power  spectrum  F(co)  of  the  underlying 

process  x(t).  Unlike  the  first  method,  it  can  be  used  to  construct 

aliases  for  any  given  power  spectrum  (or,  equivalently,  for  any  given 

correlation  function).  We  begin  by  writing  (see  [3],  p.  5?) 

27i(k+l)p 
c  (n)  -  C(nh)  »  ^^       I  exp(iconh)F(a))da> 

k=-oo  2Tikp 

(6)  /"p  ^ 

-     I       exp(i»nh)j    \        F(co+2nl<p)  L  dio     , 
-np  I  k»-oo  J 


the 


By  the  convex  hull  of  a  set  of  functions  f  ,  1  <  n  <  N,  is  meant 

N  N 

set  of  all  functions  of  the  form  \  a  f  ,  0  <  sij,  <  1>  where  \.      a^  =  1. 

n=l  n=l 
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where 

P  =  1/b 

is  the  sampling  rate.     At  this  point,   it  is   convenient  to  introduce 
the  linear  operator  y/  =Yy        ,  defined  by 

(7)  yffM   =y    f  (co  +  2nkp) 

k=-oo 

The  following  properties  of  the  operator  "Yy  are  easily  verified: 

1.  If  f(co)^L  ,   thenYyf(o))  exists  for  almost  all  co  and  has 
period  2np     . 

2,  If  f(ai)6L"'",   then>/f(a))£L"'"(-np,iip),   and  in  fact 

/IXp  /OO 


J     7/f(co)d£o  =  I       f  (co)dco 


'-np  '-co 

3.     If  f(co)  is  continuous  and  f(co)   =  0  ( Icol"-*-"®),   e  >  0,  then~yy 
f(a3)  is  continuous     • 

li«  "Yi  f(co  +  2 imp)  <=flfi(ji),  where  n  is  any  integer. 

S,     If  h(oi)  hasperiod  2iip,  then~yYh(co)f(co)   =  h(co)y/f (t^). 

Suppose  now  that  f(co)  >0,   as  in  the  case  of  power  spectra.     Then 
f(co)  can  be  regarded  as  a  distribution  of  mass  on  the  oj-axis,   anddTf (co) 
is  the  mass  distribution  on  the  circumference  of  a  circle  of  radius  p 


The  letter lOf  is  meant  to  suggest  "winding"  (see  below). 
In  this  connection,  note  that  Wiener's  class  M, ,   defined  in   [U] ,  p, 

73 >  is  essentially  the  class  of  continuous  L  functions  whose 
'•'wound-up"  versions  are  finite  (and  hence  continuous). 
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which  results  when  the  co-axis  is  wound  up  on  the  circle  in  such  a 
way  that  the  point  co  =  0  on  the  line  coincides  with  the  point  (0  =  0 
on  the  circle J  in  this  construction,  superimposed  mass  simply  adds. 
Comparing  (2),  (6)  and  (?),  we  see  at  once  that  the  relation  between 
the  power  spectrum  fj^(«^)  of  the  random  sequence  x(nh)  and  the  power 
spectrum  F(co)  of  the  underlying  process  x(t)  is  just  fj^(")  =  pHFCpco), 
i.e.,  fy^(")  is  the  "wound-up"  version  of  F(co),  measured  in  anr^lar 
units.  Thus,  any  other  power  spectrjm  F(co)  which  gives  the  same  mass 
distribution  when  wound  up  on  a  circle  of  radius  p  is  an  alias  of 
F(co),  with  respect  to  periodic  sampling  at  rate  p,  and  conversely, 
all  aliases  of  F(co)  wind  up  to  the  same  mass  distribution.  From 
this  point  of  view,  we  can  create  any  number  of  aliases  of  F(co)  by 
shifting  some  or  all  of  the  spectrum  F(co)  through  integral  multiples 

00 

of  2np.  More  precisely,  given  any  decomposition  F(co)  =  \      F  («), 

n=il 

whei^  the  F  (cti)  are  all  non-negative,   construct  the  new  spectrum 
n 


A 

where  the  k     are  arbitrary  integers.     Then  F(co)  is  even,  non-negative 
n 

and  integrable,   and  moreover 

7f"F(co)   =  I  y  |lr^^((o  +  2nk^p)  +7/F^(^  +  2Ttk^p)l 
n=l  C  J 

r 


^  V"U^V"^   +>fF^(-co)"|  -  |Wf(co)  +7/f(-co)   -VfC")   , 
n=l  I  J  L  -J 


as  required. 
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We  note  in  passing  that  the  technique  developed  above  can  also 
be  used  to  constnict  correlation  functions  which  vanish  at  all 
points  of  the  sequence  t  =  nh,  n  4  0.  For  suppose  F(co)  is  any 
power  spectrum  for  whichof  F(co)  has  a  positive  lower  bound  •  Then 
the  power  spectrum  H(co)  =  F(co)  /3TT(co)  obviously  satisfies  «H(w)  =  1, 
whence  it  follows  from  (2)  that 


c,  (n)  =  C(nh)  =  2np6 
h  "^  on 


where  6   is  the  Kronecker  delta,  i.e.  C(nh)  =  0  for  n  ^  0,  as  required, 

3.  Jittered  periodic  sampling.  When  trying  to  implement  the  periodic 

sampling  scheme,  we  have  to  expect  deviations  of  the  samp!ling  times 

from  the  valines  nh,  due  to  jitter  in  the  device  generating  the  time 

markers.  We  now  inquire  as  to  the  effect  that  such  inaccuracy  has 

on  aliasing;  in  particular,  it  is  conceivable  that  the  jitter  actually 

present,  or  some  extra  jitter  deliberately  introduced,  might  eliminate 

the  aliasing  encountered  with  ideal  periodic  sanpling.  This  leads  us 

to  consider  the  following  modification  of  ideal  periodic  sampling: 

Let  the  sampling  times  be  of  the  form  t  =  nh  +  Y„>  n  =  ...,  -2, -1,0,1, 2, ...^ 

where  the  y  sre  a  family  of  independent,  identicaU.y  distributed 

2 
Gaussian  random  variables,  with  mean  zero  and  variance  C  ,  i.e.  with 


r 

In  particular,  this  implies  that    log  fy,(<^)  dco  >  -oo,  so  that  the 
random  sequence  x(t  )  is  indeterministic  (see  e.g.  £3},  p.  69). 
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coDunon  probability  density 

The  conditj-on  for  small  jitter  is     cT  «  h;  in  the  limit  a' =  0 ,  pit) 

reduces  to  5(t),   the  Dirac  delta  function,  and  we  hnve  ideal  periodic 

sampling.     Note  that  Et     =  nh,   so  that  the  average  spacing  between 
sampling  times  is  still  h. 

Since  there  is  no  statistical  dependence  between  the  y     and  x(t), 
the  zero-mean  random  sequence  x(t   )  is  stationary''  (with  respect  to 
the  product  statistics  of  the  f     and  x(t)    )  with  correlation  function 

(8) 


.00 


/. 


cCr)  Pj^(r)  dr        , 


00 


where  p  (t)  ,  the  probability  density  of  t    -  t  ,  is  independent 
n  iTi'»"ri    in. 

of  m.  For  n  =  0,  we  have  P  (t")  =  SCt),  and  (8)  reduces  to  c,  (O)  =C(0), 

a  relation  which  is  obvious  from  the  fact  that  the  mean  square  of  any 

sequence  sampled  from  x(t)  must  be  the  same  as  the  mean  square  of  x(t) 

itself. 


See  [3],  p.  58,  The  slight  abuse  of  notation  in  (8)  is  transparent, 
ioe,,  the  first  E  is  an  average  over  the  product  statistics  of  the 
Y  and  x(t),  whereas  the  second  E  is  an  average  over  the  statistics 

of  the  Y  and  the  third  E  over  the  statistics  of  x(t). 
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Now  if  g('t')  is  the  probability  density  of  a  Gaussian  random 
variable  with  mean  zero  end  variance  2o   ,   i.e. 


(9)  g('^)  =    -i exp  (-t^Ao^) 

2  -/rTcr 


we  see  at  once  that 

r 

(ICO  c^(n)  -   I       g(nh-r)  C(t)d'r    ,     n  +  0 

Since  g(nh  -t  )  does  not  reduce  to  6('C)  for  n  =  C,(10)  is  not  valid 
for  n  =  Oj   instead,  as  already  noted,   we  have  c,  (0)   =  C(0).     However, 


we  shall 

assume  that 

the 

quantity 

,00 

(11) 

d,(0) 

1 

-00 

giv)C{t)dt, 

obtained  by  setting  n  =  0  in  the  right-hand  side  of  (10),  is  known 
fron:  a  suitable  side  experiment.  For  example,  we  might  sample  not 
only  at  the  points  t  ,  but  also  at  the  points  t     +  h,  and  calculate 

(12)  Vo^-.T^^    iV"    ^K^^^^K^l^      ' 

"  N->oo  N  /  ■  n  n+1 

Or  we  might  mark  off  two  sets  of  random  tines  t  =  nh  +  v  , 

n       'n' 

t  =  nh  +  Y^  >  where  the  Y„  are  another  set  of  independent  Gaussian 


There  is  perhaps  a  slight  contradiction  in  (12),  for  although  the 
times  t  are  subject  to  jitter,  it  is  assumed  that  no  additional 
jitter  occurs  in  marking  off  the  times  t  +  h. 
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random  variables  with  common  distribution  (9)>  and  calculate 

n=l 


(The  c.  (n),  n  4  0,   can  be  calculated  from  either  the  t     or  the  t   ,) 
h  '  n  n 

Our  primary  justification  for  assuming  that  d.  (O)  is  known  is  that 
it  is  easily  calculable  from  an  ensemble  point  of  view.     In  fact, 
indicating  a    "  tj-pical"  infinite  set  of  sample  functions  of  the  process 
x(t)  by  x.(t),   1  <  i  <  00,  we  have 

The  issue  involved  is  that  the  various  points  labelled  t     (n  fixed) 
are  scattered  in  time  when  looked  at  across  the  ensemble.     In  any  event, 
since,  as  we  shall  see  below,  aliasing  is  in  general  present  even 
when  d,  (O)  is  known,  it  is  present  a  fortiori  in  the  absence  of  this 
extra  information. 

We  observe  next  that  the  convolution  of  C(t)  with  git),   i.e. 

,00 

(13)  D(T)   =      1        g('C-u)C(u)du 

-00 

is  itself  a  correlation  function,  as  follows  st  once  from  the  fact 
that  it  has  the  non-negative  Fourier  transform  H(a))  =  F(co)G(co),  where 


G(w)  =  exp(-o^ci)  /2) 
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Thus  the  sequence 


\M     =     c^(n)        ,  n  +  0       , 


\{0)     =     J        g('t)C(T)dt'  , 

/Loo 

is  clearly  a  correlation  sequence.     In  fact,   if  y(t)  denotes  a 

zero-mean  random  process  with  correlation  function  Biv),  then  d,  (n) 

is  the  correlation  function  of  the  random  sequence  y(nh)  obtained 

by  periodic  sampling  of  y(t)  with  spacing  h.     Incidentally,   this  gives 

an  independent  proof  that  the  sequence  c,  (n)  is  a  correlation  sequence, 

since 

c.  (n)     -  d.  (n)  +  a6 
h  n  on  , 


where 


r 

]       F((»)  {i  -  G(co)y   dco  >  0 
/Loo 


Thus,  c,  (n)  is  the  sum  of  two  correlation  functions  and  therefore 
itself  a  correlation  function. 

The  question  of  whether  or  not  aliasing  occurs  with  Jittered 
periodic  sampling  can  now  be  stated  simply  as  follows!     fire  there  two 
(or  more)   correlation  functions  which  take  the  same  value  at  't  =  0 
and  which  when  convoled  with  g(T)   give  the  same  values  on  the  lattice 
'^  =»  nh,   n  »...,  -2,-1,0,1,2,...?    Or,   in  spectral  terms,  are  there 
two   (or  more)  power  spectra  which  have  the  same  total  mass  and  which 
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when  multiplied  by  G(co)   give  the  same   "wound-up"  version  on  a  circle 
of  radius  p?     (Note  that  both  statements  reduce  to  the  case  of  periodic 
sampling  li^en  cf  =  0.)     As  we  shall  now  show,  the  answer  is  that  aliasing 
persists  in  the  presence  of  jitter  (with  an  exception  to  be  discussed). 
We  now  proceed  to  give  two  ways  of  constructing  aliases,  which  parallel 
the  methods  used  before  for  constructing  aliases  in  the  case  of  ideal 
periodic  sampling. 

The  first  method  consists  of  finding  two  correlation  functions 
D  (f)  and  D  (T)       which  vanish  on  the  lattice  't  =  nh,   n  =...,-2,-1,0,1,2,. 
and  whose  Fourier  transforms  have  the  same  integral  when  divided  by 
G(co).     Since  G(aj)  ->  0   as  co  ->  +oo,  we  shall  arrange  for  H_j_(co)  and  H_((b), 
the  Fourier  transforms  of  B^^it)  and  D_(ir),  to  vanish  outside  of  some 
finite  interval.     This  can  be  achieved  by  choosing  C_j.(r)  and  D_(tr) 
to  be  entire  transcendental  functions  of  exponential  type,  whence  it 
follows  by  a  theorem  of  Paley  and  Wiener  (see  e.g.    [^   p.  13U)  that 
H  (co)   and  H  (oo)  vanish  outside  of  some  finite  interval. 

With  this  in  mind,  define  the  function 

1-   f^\      I  sin  Tcf  1        .   2     ^ 
\{V)  =1 I     sm  npf  , 

where  a  is  an  arbitrary  real  number  and  p  =  lAi  is  the  sampling  rate. 
Then  form  the  linear  combination 


h('Zr)  =  A  h^C-C)  +  hp(0    , 


The  notation  anticipates  the  nature  of  the  construction. 
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where  a  4  P  ar^d  A  is  a  constant  to  be  chosen  later.  The  function  h(r) 

1       2 

is  rea] ,   even  and  in  L  /I  L  ,   and  vanishes  at  all  points  of  the  lattice 

f  =  nh.     Its  Fourier  transform  H(aj)  is  therefore  real,   even  and  in  L  . 
Moreover,  h(T)  is  an  entire  function  of  exponential  type  with  exponent 
CO     *  2n(p  +  max(a,p)).     Hence,  by  the  Paley-V7iener  theorem,   H(ai) 
vanishes  for   lw|  >  co  ,   so   that  both  H(cci)  and  H(a))/G(a))  are  in  L  , 
We  now  choose  the  constant  A  so  that 

,00 

H(a) 


GfcoJ  ^  =  ^ 


/-OD 


Since  the  functions  h  (f )  and  hpCt')  are  linearly  independent,  the 
resulting  h(f)  cannot  vanish  identically.  Now  define  the  non-negative 
functions  H  (oo)  and  H_(cii)  exactly  as  in  equations  (3)  and  (U)  of  the 
preceding  section,  and  let  their  Fourier  transforms  be  D  (f )  and 
C  (f  )>  which  are  automatically  correlation  functions.  Then  D  (f )  and 
D  (t)  agree  on  the  lattice  "V  =   nh,  and  moreover,  by  construction  the 
corresponding  correlation  functions  C  (t)  and  C  (f),  defined  by  (13), 
agree  for'C=  0.  This  completes  the  construction  of  aliases  by  the  first 
method. 

The  second  method  of  constructing  aliases  is  based  as  before  on 
the  winding  construction,  and  has  the  advantage  that  it  allows  one  to 
produce  aliases  for  any  given  power  spectrum.  Using  the  correlation 
sequence  d,  (n)  we  form  the  function 

00 


H(co)  =  ^  y-'  exp(-ia)nh)dj^(n)  =  F(co)G(q)  ,  oa^N  , 
n=  -00 

H(co)  =  0  ,  ciiN, 
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where  N  is  the  Nyquist  band,  defined  by   jcoj   <  np.     We  also  have 
available  the  quantity  c,  (O )  =   |     F(a))dci>,     Now  evaluate  the  integral 


P 

-00 


T         I  H(6j)    , 


We  have 


^T^J  ST'  F(u+  2npn)G(co  +  2npn)  l  dm 


n=-oo 


-j      J  21    F(co  +  2npn) 
'^Tip  n=-oo 


fCO 


dt£>     =      I       F(a))dco  =  c,  (C) 


» 


'-00 


where  the  equality  holds  if  and  only  if  F(co)  lies  entirely  in  N  (i.e. 
vanishes  almost   everywhere  outside  of  N).     Thus,  if  the  values  of 
I  and  c.  (O)  are  equal,  we  know  that  F(co)  lies  in  Nj   consequently, 
F(co)  =  H(co)/G(co)  and  we  can  determine  F(co)  uniquely  in  this  case. 
Moreover,  unlike  the  case  of  ideal  (i.e.   jitter-free)  periodic  sampling, 
we  can  ascertain  whether  or  not  FCu)  lies  in  N  by  a  simple  test,  namely 
comparison  of  I  and  c.  (C),   instead  of  having  to  know  a  priori  that 
F(co)  lies  in  N. 

If  F(w)  does  not  lie  in  N,   then  I  <  Cy-(0)»  and  aliases  of  F(co) 
can  be  constructed  as  follows.     Define  the   "virtual  spectrvun" 

F(»)  =  H(co)/t}(co)       ,  co€N     , 

A 


F(o))  =0  ,  co^N     J 
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F(co)  =  F(6o)   if  and  only  if  T{oi)  lies  in  N.     Now  write  F(co)   in  any  way 
as  the  sum  of  two  non-negative  functions  F- (co)  and  Fp(6)),  and  define 

;pj— r  ■j  F-(co  +  2itpn)G(u  ♦  2npn)  ♦F^Cna  +  2npn)G(-a  +  2npn)|  , 


*   2G( 


where  0  <  a  <  1  and  n  is  arbitrary,      (Note  that  F       (u)  =  F(co),)     Then 
—      —  n,o 

F       (to)  is  even  and  non-negative,  and  moreover 


n,a 


V/f„  „(co)G(co)   =  7fF(co)G(co)   =  H(co)        . 
Since  Fp(aj)  vanishes  for  co  C  N,   it  is  easily  seen  that  if  n  4  0 

^7— T    i'^2^^  ♦2npn)G((i5  +  2npn)  +F2(-co  ♦2npn)G(-co  +  2npn)>'  dw 


/OO 

>-00 


>    1    J  F^Cw)  +  ^g^^^y  '^     » 


i.e.,  the  distribution  F   (co)  has  more  mass  than  the  initial  distri- 

n,a 

bution  F(co)  if  0  <  a  <  1.  In  fact,  since  the  first  integral  increases 
with  n,  there  exists  an  integer  n  >  0  such  that 

^00  /OO 

F     ,(ai)  dw     >  F(co)dco 

/-  00  '-00 


■«•, 


A 


We  exclude  the  trivial  case  F,(co)  =  0  or  Fp(co)  s  0, 
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for  all     n  >  n  .     Therefore,  since 


00  rOO  j<X) 


F(co)dco  <        F(oi)dco     , 


F       (a3)do3  = 
n,o  J  J 

'-00  ''-00  -  00 

and  since  F       (co)  depends  continuously  on  the  parameter  a,  there  is  a 
njO 

value  of  a,   say  a   ,  between  0  and  1,   such  that 

/CO  /OO 

F         (co)da5     =      I     F(cLi)dco  . 

;     "'%  )  ^ 

-00  -00 

It  is  clear  that  any  such  function  F    _    (o))   is  an  alias  of  F(aj)  with 

'  n 

respect  to  ^iittered  periodic  sampling.  This  completes  the  construction 

of  aliases  by  the  second  method. 

li.  Additive  random  sampling. 

The  reason  that  aliasing  is  still  present  with  jittered  periodic 
sampling  seems  to  be  that  the  sampling  tines  t  are  still  "  attracted" 
to  the  equi-spaced  values  nh.  What  seems  to  be  needed  to  break  up  this 
regularity  is  some  sort  of  "floating  point"  sampling  scheme.  Thus,  we 
consider  next  additive  random'  sampling,  in  which  each  sampling  time  is 
derived  from  the  preceding  one  by  the  addition  of  an  independent  random 
variable*^.  Specifically,  let  t^  =  t^.^+Y^j  »  where  the 
Y  ,  n  «  ...,-2,-1,0,1,2,..,,  are  a  family  of  identically  distributed. 


The  considerations  of  this  section  can  be  extended  to  the  case 
where  gCtr)  instead  of  being  Gaussian  is  a  function  whose  Fourier 
transform  has  isolated  zeros  only,  e.g.  when  g(T)  vanishes  outside 
of  a  finite  interval.  Again  it  is  found  that  aliasing  persists.  It 
is  also  clecr  how  to  extend  the  second  method  of  constnicting  aliases 
to  give  aliases  of  a  more  general  type, 

If/ith  jittered  periodic  sampling,  the  differences  between  successive 
sampling  times  are  not  independent  random  variables,  and  the  variance  of 
t  is  20*2  for  all  n.  On  the  other  hand,  with  additive  random  sampling, 
the  variance  of  t  goes  to  oo  with  n,  if  sampling  starts  in  the  finite 
past  (whence  the  phrase  "  floating  point"  ),  although  Et^  still  equals  nh. 
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independent  random  variables,  with  Ey„  =  h  <oo  and  comnon  probability 

"  ^00 

density  pCx:).  Of  course,  pCir)  >  C  and  I  p(r)d'C' =  1,  \-Je   assume  that 

2  -00 

p(f)  is  also  in  L  and  that  p{f)   =  0  for  'Zr'<  0,  The  last  condition 

corresponds  to  the  reasonable  requirement  that  a  sample  time  with  a 

given  index  should  come  after  all  sample  times  with  smaller  indices. 

For  later  use,  we  introduce  the  Fourier  transform 

^00  /OO 

(lb)    ^(co)  =  J  exp(irot)p(tr)dr  =  J  exp(i^f)p(tr)dtr    , 
-00  '^o 

the   characteristic  function  (in  the  sense  of  probability  theory)  of 
the  distribution  p{v) ,     The  function  ^(co)  is  uniformly  continuous  on 
the  whole  real  line  and  satisfies  the  conditions 


0(0)   =  1   ,      |0(co)|   <  1  ,  (/(^)   =  WU^,    ^^^^  0(co)   =  0      , 

where  the  overbar  denotes  the  complex  conjugate.  It  will  be  important 
to  extend  (lU)  to  complex  arguments  by  writing 


,00 


(15)  iZf(s) 


exp(ist)p(f)dr 


/. 


where  3=0  +  io)'.     Eq.(l5)  defines  a  function  which  is  analytic  and 
bounded  in  the  open  half -plane  Im  s   >  0,   and  continuous  in  the  closed 
half-plane  Im  s    >  0.     Moreover,   for  Ims>  0,  we  have   |?'(s)|   <  1,  with 
equality  holding  only  for  s  =  0, 

As  in  the  case  of  jittered  random  sampling,  the  correlation  function 
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of  the  random  sequence  x(t  )  is  given  by 


c,  (n)  =  Ex(t     „)x(t   )   =  E  Ex(t     „)x(t   )i  t     „  and  t^  fixed 
n  m+n        m  I         m+n         m        m+n  n  j 


(8) 

/OO 


c(f)Pj^(tr)d'rr    , 


where  now,   in  view  of  the  new  definition  of  t  ,  the  prob3bility  density 
p  (t')  is   given  by 

r  C 

PjjCtr)  =  p^_^(f-u)p(u)du       =     I   p^_^(r-u)p(u)du     ,     n  >  2  , 

(16)  <-oo     "  -^       ' 

P^(tr)  =     pCf)     , 

i.e.,  p  (f ),  n  >  2,   is  obtained  by  successively  convolving  p(f )  with 
itself  n  times.     As  before,  p  (XT)   =  6(f)  and  c.  (-n)  =  c,  (n).     Note 
that  p  (T),   n  >  1,  vanishes  for  negative  f  ,   so  that   (8)  becomes 

~/00 


c^(n)  = 
(17)  ^ 


C(t)p^('6-)df      ,  n  >  1     , 


Cj^(O)  -  C(0)        . 


Moreover,  using  Parseval's  theorem,  we  have 


^00  /OO  /OO 

2,^N.^  I      \rj,   M2n,     ^     I      \rif.\\2^  I  _2 


p^('K)dt    =     I      |9^(co)|^"dfo  <    j      |9f(co)|'^da)  =     I   p'^('?r)d^      , 
-  00  -  00  'o 


2 
so  that  p  (tr),   n  >  1,  is  in  the  class  L  .     Again  applying  Parseval's 

theorem  and  the  convolution  theorem  for  Fourier  transforms,  we  find 
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that  the  spectral  equivalent  of  (1?)   is 


fOO 


(18)         Cj^(n)     =  F(co)9f"(co)  do)     ,         n  >  0 


-00 


The  reality  of  c,  (n),   as  defined  by  (18),   follows  from  F(-co)   =  F(co) 
and  <^i-ci)   =  j5Tco7. 

Equations   (17)   and  (18)  exhibit  a  direct  relation  between  the 

correlation  sequence  c,  (n)  and  the  correlation  function  C(t>)  or  power 

gsectrum  F((a),   and  we  shall  use  these  formulas,   particularly  (18),  to 

study  aliasing  in  the  case  of  additive  random  sampling.     It  is  of  some 

interest,   however,   to  examine  the  relation  between  the  spectrum  F(w) 

and  the  spectinm  f^(w)  of  the  random  sequence  x(t   ).     To  do  this,  we 

use  (2)   and  (18)  to  write 

/OO  ,  _ 

f^(<o)=-|^^    exp(-icon)   c^(n)   =|5jJ     F(u)|i +^  0"(u)(e^"^-.e-Wdu, 
n=-oo  -00  n=l 

where  the  second  series  is  convergent  for  all  u  I'  0 .     Carrying  out  the 
suitmstion  and  using  the  fact  that  F(a))   is  even,  we  obtain 


fh(->  ■  h 


(1?) 


„  ll-e    0(u), 


.  1  Cnu)  llzMs)£ 

2«J.„  ||l-e-0(u 


so  that  fr^('^)  exists  and  is  non-negative  everywhere  in  (-n, it),   except 
perhaps  at  co  =  0,  where   (19)  becomes  indeterminate.     Incidentally,   this 

*See   [3],  p.  58. 
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gives  an  independent  proof  of  the  fact  that  c   (n)  is  a  correlation 
sequence,   i.e.  that  the  sequence  x(t   )  is  staticnaiy.     Finally,  we 


have 


,n  /OO  f  /IT 


/-n 


f^(co)  dco  =  ^   1     F(u). 


fl^e^(u)l 
'll.3^i2f(u)/ 


Rer       i.!^^      )  duy  du 
»•  -n 


-i-         F(n)JRe  ^  *  '''<"'     Mdu 

\^  |z|=l  / 

as  is  to  be  expected. 

The  question  of  whether  aliasing  occurs  with  additive  random  sampling 

can  now  be  simply  stated  as  follows:     For  which  p(t)   is  there  only  one 

correlation  function  which  leads  via  (18)  to  a  given  correlation  sequence 

c,  (n)   ?    Or  in  spectral  terms,   for  which  p(r)  is  there  only  one  power 
n 

spectrum  which  leads  via  (19)  to  c,  (n)  ?  We  now  prove  a  theorem  which 
gives  a  partial  answer  to  this  question,  i.e.  which  gives  a  sufficient 
condition  for  additive  random  sampling  to  be  alias-free. 

Theorem  1.     Additive  random  sampling  is  alias-free  if  the  characteristic 
function  j6(co)   takes  no  value  more  than  once  on  the  real  axis. 

Proof.     If  aliasing  is  possible,  then  (18)  is  satisfied  for  two 

1   2 

distinct  real  ( no n- negative)  functions  in  L  f\L   ,  i.e.,  there  exists 

12  * 

a  i^al  function  H(co)  in  L  Pi  L  >  which  is  not  a  null  function  ,   such 

that 

(20)  J       H(co)0"(co)dco     =     0 

^-00 


By  a  null  function  is  meant  a  function  which  vanishes  except  on  a 
set  of  measure  zero. 
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for  n  >  0.     Thus,  we  must  show  that  if  J?5(a))  is  one-to-one  on  the  real 
axis,   then  any  H((X))  satisfying  (20)  must  be  a  null  function. 

We  begin  by  showing  that  if  (/(co)  is  one-to-one  on  the  real  axis, 

then  0{u>)  is  actually  one-to-one  for  Ims  >  0,  or,   in  the  terminology 

of  conformal  mapping,  0(s)  is  schlicht  in  the  closed  upper  half  plane. 

To   see  this,  we  note  first  that   it  follows  by  a   Phragra^n-Lindelof  theorem 

(see    [6],   p.  179)   that  the  boundedness  of  0(a))  together  with  0(co)  ->  0 

as  03  ->  +00  implies  that  ^(s)  ->  0  uniformly  in  the  upper  half-plane, 

i.e.   that  0(s)   can  be  made  as  small  as  we  please  outside  a  suitably 

large  semicircle   P^   erected  on  some  diameter  (-P.,P).     Let   P  be  the 

closed  contour  consisting  of   p,   and  the  segment   (-R,R).     Now  if 

a  ^  0  is  a  value  taken  by  0(s)  at  some  point  in  Ims  >  0,   choose  P  large 

enough  so  that    |S2^(s)l<a  on  P.  By  the  argument  principle   (see   [S] ,  p.  116), 

if  0(s)   takes  the  value  a  more  than  once  for  Ims  >  0,   then  the  image  of 

the  contour   P  under  the  mapping  z  =  0(s)  would  have  to  wind  more  than 

once  around  the  point  a.     However,   this  is  impossible,   since  by  hypothesis 

the  image  of  the  segment  (-R,R)  of  the  real  axis  cannot  intersect  itself, 

whereas  the  image  of   p^    is  confined  to  a  circle  of  radius  less  than  a. 

Moreover,   if  0(s)  takes  the  value  a  at  some  point  s     of  Ims  >  0,   it 

cannot  take  the  same  value  at  a  point  co     of  the  real  axis,   for  since 

o 

0(s)  takes  all  values  sufficiently  close  to  a  in  a  neighborhood  of  s 
(see    [7],  p.  I8U),  and  since  0(s)  is  continuous  at  co  ,  we  could  then  find 
two  non-zero  values  near  a  taken  twice  in  the  open  upper  half-plane,  which 
is  impossible,   as  we  have  just  seen.     It  follows  that  (^{s)  can  take  no 
non-zero  value  more  than  once  in  Ims  >  0 . 
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Moreover,  )6(s)   cannot  vanish  for  Iras  >  0,   for  if  ^(s   )  =  0,   Ims     >  0, 
then  in  a  neighborhood  of  s     which  lies  entirely  in  the  open  upper 
half -plane  j?y(3)taices  all  sufficiently  small  values.     But  since  0(s)->  0 
as  s  ->oo,  we  could  then  'find  a  small  non-zero  value  taken  twice  in  the 
open  upper  half -plane,  ^ich  is  impossible,  as  already  noted.     Moreover, 
0(a))   cannot  vanish  for  real  co,   since  then  ^(-co)  also  vanishes,   which 
violates  the  hypothesis  of  the  theorem  (recall  that  0(0)  =|=  0.)     Thus,  we 
have  finally  shown  that  0(s)  is  schlicht  in  the  closed  upper  half  plane. 

We  now  turn  to   the  core  of  the  proof.     Consider  the  analj-iiic 
function  z  =  /5(s),  which  maps  the  half-plane  Ims  >  0  onto  a  domain  J /, 
in  the  z-plane  which  is  entirely  contained  in  the  circle    lz|   =1  and  whose 
boundary  does  not  intersect  itself  and  passes  through  the  points  z  =  0 
and  z  =  1.     Next  consider  the  analytic  function 

^=4(s)=        ^ 


s  -  B  +  i 


where  p   is  a  real  constant}  C(s)  maps  the  half -plane  Ims  >  0  onto  a 
circlel  Ly,  in  the  ^-plane  which  has  5^  =  0  on  its  boundary.     Now 


X.    =    lis)    =    l{r\z))    H    J(Z) 


intl  L 


is  a  single-valued  fvmction,  analytic  in  the  domain a  L    and  continuous 
on  the  boundary  of  J  L  ,  which  mapsikij     ontoAiy  in  such  a  way  that  the 
points   z  =  0  and  X.  =  0   coincide.     By  a   theorem  of  Walsh  (see    [S]  ,   p.   36), 
given  any  e  >  0,   there  is  a  polynomial  P(z)  =  c^  +  c^z  +   ...  +  c^z 


*By  0~   (z)  we  mean  the  inverse  function  of  0(z). 
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such  that 


|j(z)  -  P(z)|  <  c 

on  the  boundary  B  of i  L  .  Since  J(0)  =  0,  we  have  |P(0)|  <  e.  Thus, 
setting  P,(z)  =  P(z)  -  c^,  so  that  P^(0)  =  0,  we  have 

|j(z)  -  P^(z)|  <  2e 


on  B   .     It  follows  that 
z 


(21) 


U(co)  -  P^(0(co))|   <  2e     , 


since  0(s)  maps  the  real  axis  into  B  * 

Suppose  now  that  H(co)  satisfies  (20 )  for  n  >  1,  Then,  using  (20) 
and  (21),  we  have 


00 


H(co) 


/ 


CO  -  p  +  i 


r  dw 


00 


<00 


H(w) 


-00 


|4(co)-P^(^(co 


m  cao) 


^00 


<  2e  I   |H(co)|dfo. 

-CD 


Since  K(co)CL     and   e  is  arbitrary,   it  follows  that 

rOO 


(22) 


H(co) 


■'-co 


CO  -  p  +  i 


d»  =  0     . 


Taking  the  insginarj'-  part  of  (22),  we  obtain 

/OO 


H(co) 


/- 


00 


1  +   (co  -  p)' 


dco     =     0     , 


since  H(co)  is  real.  Recalling  that  p  is  an  arbitrary  real  number,  we  see 
that  the  convolution  of  H(co)  and  l/(l  +  a>  )  vanishes  identically.  It 
follows  from  the  convolution  theorem  for  Fourier  transforms  that 
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h(f  )exp(-|T'| )   is  a  null  function,   where  hCf)  is  the  Fourier  transform 
of  H(co),     But  since  expC-l^D  does  not  vanish  for  finite  f,  hiv)  itself 
and  finally  H(ai)   is  a  null  function,   i.e.   aliasing  is  impossible, 
(Note  that  because  of  the   special  nature  of  the  mapping  J(z),  we  have 
not  used   (20)  for  the  value  n  =  0,   so  that  knowledge  of  the  number 
c,  (0)   is  superfluous.)     This  conDletes  the  proof  of  Theorem  1, 

Next  we  prove  a  conditional  converse  of  Theorem  1. 

Theorem  2.     If  the   characteristic  function  0(s)   takes  the  same  value 
at  two  different  points  of  the  open  upper  half-plane,  then  aliasing  occurs 
with  additive  random  sampling. 

Proof.     We  have  0(s,)  =  ^(Sp)*  where  s^  k  Sg  and  Ims^  >  0,   Ims^  >  '^. 
It  is  asserted  that  there  exists  a  non-null  real  even  function  in  hflL 
such  that  (20)  is  valid  for  n>0.     Equivalently,       if  h(r)  is  the 
Fourier  transform  of  H(oj)  and  fj  is  the  class  of  Fourier  transforms  of  L 
functions,   it  is  asserted  that  there  exists  a  non-null  real  even  function 
h(f)  inflriL^,  with  h(0)  = 


/-co 

H(cii)dco  =  0,   such  that 


.         ^-00 
rOO  /OO 


(23)  I     h('t)p^('C)dt     =     I  h(r)p^('^)dr    =       0     ,       n  >  1     , 


'-00 


where  the  p  iv)  are  given  by  (l6).     We  shall  show  that  the  assumption  that 
there  are  no  such  H(co)  and  h(f )  leads  to  a  contradiction,  thereby  proving 
the  theorem. 

Suppose  there  is  no  h(t)  with  the  specified  properties  which 


Recall  that  Fourier  transformation  is  a  unitary  transformation  of  L     onto 
itself  (see   [U] ,   p.  70). 
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satisfies  (23).  Then  since  the  class  /=^  of  such  h(t)  is  dense 
in  the  real  function  space  L  (0,oo),  it  follows  that  there  is  no 


no 


n-null  function  in  L  (0,oo)   for  which  (23)   is   satisfied,   and 


'J 

consequently  that  any  fxmction  in  L  (0,od)  can  be  approximated  in 
L  -norm  by  linear  combinations  of  the  Pj.('^)j  "^  ^  (see  [.^ J  >  P»  37). 

Equivalently  (see  footnote,  p.  25),  the  Fourier  transform  of  any 

2  2 

function  in  L  (0,oo)   can  be  approximated  in  L  -nDrm  by  linear 

combinations  of  the  powers  4  (w)^   n  >  1.     We  choose  this  function  in 
L  (0,oo)  to  be  exp(-'J'),  with  Fourier  transform  l/(l-iio)j   this  choice 
is  motivated  by  the  fact  that  l/(l  -  is)  is  a  schlicht  function.     Then 
we  liave 


(2U) 


lira 
n->oo 


/OO 


/ 


00 


P  (;^(co))   -        ':^ 


dco     =    0  , 


where 


2  n 

P  (x)  =  a  ,x  +  a  „x     +   ...  +  a x 

n  nl  n2  nn 


is  a  suitable  polynomial  with  real  coefficients. 


Now  form  the  integral 


(25) 


1 

753: 


T7-{V^(^))-r^}^^ 


r 


To  see  this,   note  that  the  class  of  step-functions  is  dense  in 

L  (0,oo)   (see   [h] ,  p.   2U)  and  that  each  step  can  be  approximated  by  a 
"-trapezoidal  function"  ;   in  particular,   a  step  at  the  origin  can  be 
approximated  by  a  trapezoidal  function  vanishing  at  the  origin. 

Since  every  trapezoidal  function  is  inflO  L     (see   [h J ,  p.   89),  a 

subset  of  )k^  is  dense  in  L  (0,oo).      (Here,  and  below,  we  depart  from 

our  usual  notation  and  use  L  (0,oo)  to  denote  the  space  of  real  ftinctions 

h(tr)  for  which        h^(i:r)df  <  00.    ) 
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where  Ira  s     >  0  and   P  is  a  contour  like  the  one  used  in  Theorem  1, 
o  ' 

i.e.  a  large  semicircle  P  erected  on  the  segment  (-R,R)  as  diameter. 
By  Cauchy's  theorem,  the  integral  (25)  is  just  the  expression  in  curly 
brackets  evaluated  at  s  =  s  .  Furthermore,  it  is  easy  to  see  that  the 
contribution  to  (2$)  coming  from  P,j  the  semi-circular  part  of  the 
contour  P  ,  goes  to  zero  as  R  ->oo,*  Combining  these  facts  and  using 
the  Schwarz  inequality,  we  have 

/OO 


V«'(=o»-T^ 


1 

2Tii      I       CO 

'-00 


■r^lV^^'^^^-T^)^ 


(26)     <^ 


P  (i2f(co))  -  -1^ 
n  1-  iii 


•'-00 


P^(0(^))  -^ 


2a> 


where  k  can  be  chosen  so  that  the  inequality  holds  uniformly  for  s^ 
in  any  closed  bounded  subset  of  the  open  half-plane  Ims  >  0.     Choose 
this  subset  to  contain  the  points  s,   and  s     for  which  0(3, )  =  ^^(32). 
Then  by  (2h),   if  n  is  large  enough,  we  have  both 


and 


whence 


<  k  e 


<  k  e 


1  -is,   "  1  -iSp 


<  2  k6 


*  As  already  shown,  as  s  ->oo  ,  ^(s)  ->  0  uruformly  in  the  upper 
half-plane. 
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or,  since  e  is  arbitrary 


(27)  -^-X—     o     .,    ^  . 

1  -  IS,  1  -    ISp 


However,    (2?)   is  impossible  for  s-   "f  Sp   .     This  contradiction 
establishes  the  existence  of  a  real  even  H(co)  ±n  L  f\L     satisfying 
(20)   for  n  >  0,     The  constriiction  of  aliases  now  proceeds  in  the 
usual  way,   i.e.  we  introduce  the   functions  H  (co)  and  H  (co)  defined 
by  (3)  and  (U).     Then  H_j_(co)  and  H  (co),  and  more  generally  the  functions 
H  (co)  defined  by  {$) ,  are  a  set  of  alias  power  spectra.     This  completes 
the  proof  of  Theorem  2. 

From  an  abstract  point  of  view,   the  question  of  whether  or  not 
aliasing  is  present  with  a  given  method  of  sampling  is  just  the 
question  of  whether  or  not  the  mapping  from  a  space  whose  elements 
are  zero-mean  Gaussian  random  processes  to  another  space  whose  elements 
are  zero-mean  Gaussian  random  sequences  is  a  one-to-one  mapping. 
Thus,  the  fact  that  in  the  case  of  additive  random  sampling,  a  decisive 
role  is  played  by  the  question  of  whether  or  not  the  mapping 

,00 

z  =  0(s)     =  exp(ist)p(r)  dr 

is  one-to-one  is,  at  least  retrospectively,   not  too   startling. 

The  method  for  explicitly  reconstructing  the  correlation  function 
C(f)  of  the  underlying  sampled  process  x(t),   in  a  case  where  aliasing 
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is  absent  (as  revealed  by  Theorem  l)   is  straightforward.     First 
we  note  that  the  functions  ^  (co),  n  >  1,  and  correspondingly  the 
functions  p  (f ),  n  >  1,  are  linearly  independent?   for  if  a  linear 
combination  of  the  0   (co)  with  non-zero  coefficients  vanishes  identically, 
e.g«  a  .5((co)+a  ^  (©)+,, .+a     ^  (co)  a  0     then,   since  0(co)  is  continuous, 
factoring  this  polynomial  leads  to  the  contradictory  conclusion  that 
0(co)  =  const.     Next,  we  apply  the  orthogonalization  procedure  to  the 
fxinctions  p  (V),   n  >  1,  obtaining  the  orthonormal  set  q   (f),   n  >  1, 
which  obey  the  relations 

q„(^')q„('lr)df    =     6^„ 
in        n  ran 

'o 

We  then  express  the  q  (t)  in  the  usual  fashion  (see  e.g.    [9j ,  p. 22) 
in  terms  of  linear  combinations  of  the  p  (tir),   i.e. 

(28)  q^(Tr)  =  bj^PTL(^)  +  ^02^2^*^^  *   •**  *  ^nnV'^^   ,     n  >  1     . 

Multiplying  (28)  by  C(f )  and  integrating  from  0  to  co,  we  obtain  the 
formula 

(29)  Pj^  =     q^('C-)C(f)df  =  b^c^(l)  +  b^c^(2)+  ...  +  b^c^(n)>  n  >  1, 


relating  the  generalized  Fourier  coefficients  p     to  the  correlation 
sequence  c,  (n).     Finally,   in  terms  of  the  (?     and  the  q_('^)*  we  have 


00 

(30)  CCf)   =  V"  pj^qj^('C)    ,  r>0     , 

n=l 

where  the  expansion  is  valid  in  the  sense  of  mean  convergence  (i.e. 
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2  \  ^ 

conversence  in  the  L  -norm).       Eqiiivalently,  we  can  develop  an 

expression  analogous  to   (30)  for  F(co),   the  power  spectrum  of  x(t), 
necalling  that  C(-f)  =  gCt)   and  taking  the  Fourier  trsnsform  of  (30), 
xre  have 

/OO  /OO 

F(oi)  = -^^  C('^)exp(-iw^')df  =  l^  ^  C(f  )|exp(iu«-)+exp(-iu>'i')|df 

^00  n=l    r> 


(31) 


•J-  Y~  B  -(t  («)   +  Y  (^) 
2n  /     '^n^  n  n 


where 

(32)     Y^(w)  =  J   q(tS-)exp(ico't)df    =  b^g^(co)  +  b^^^^Cw)*..  .-^b^i^^Cco),   n  >  1 

and  the  coefficients  6     and  b       are  the  same  as  in  (28)  and   (29).     Since 

n  nm 

^(•<o)  =  0(coj   ,  we  have  1  (-co)  =  Y~(coJ,   30  that  the  even  function 
defined  by  the  last  term  of  (31)  is  real,   as  required.     Again  the 
expansion  (3I)  is  valid  in  the  sense  of  mean  convergence.     (Actually, 
the  validity  of  (30)  and  (3I)  in  the  sense  of  mean  convergence  requires 
the  completeness  of  the  functions  p  (li*),   n  >  1,   in  L  (C,oo).     This 
does  not  quite  follow  from  our  proof  of  Theorem  1,   since  there  V7e 

showed  that  (20)  implies  that  H(aj)  is  a  null-function  with  the  aid  of 

1  2 

the  assjmiption  that  H(a))  is  in  L     as  well  as  in  L  .     By  a  more   involved 

argument,  one  can  establish  the  required  completeness  of  the  p  (V)  in 
L  (0,00),  i.e.  one  can  drop  the  hypothesis  H(co)£l  .     (See    [loj). 


Since  C('6')  is  continuous  by  assumption  and  in  practice  will  generally 
have  further  regularity  properties,  i-re  can  expect  that  the  expansion 
(30)  will  usually  be  convergent,  or  at  least  suramable  by  arithmetic 
means.     Then,  since  the  number  c,  (O)  has  not  been  used  in  forming  (30 )j 

we  can  check   (30)  by  calculating  C(0)  and  seeing  whether  c.  (O)  =  C(C) 

holds.     Vfe  can  also  expect  that  (3I)  will  usually  be  convergent  or 
summable. 
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$»   Alias-free  sampling  methods. 

Using  the  theory  of  the  preceding  section,  we  now  exhibit  a 
number  of  examples  of  additive  random  sampling  methods,   some  of  which 
are  alias-free  and  some  of  which  are  not.     Consider  first   "  Pbisson 
sampling"    corresponding  to  the  choice 

p(f)  =  p  exp(-pt)   ,  t>0, 

(33) 

pCi-)  =0  ,        r<o  , 

which  by  (l6)  leads  to 


p^Ct-)  =  0  ,  't<0   , 

for  n  >  Ij  with  the  choice  (33)  the  sampling  times  are  the  occurrence 
times  of  the  events  in  a  Poisson  process  with  average  rate  p.  Since 
the  corresponding  characteristic  function 

0(s) 


P  -IS 

takes  no  value  more  than  once  on  the  real  axis   (or  anywhere  in  the 
complex  plane),  it  follows  by  Theorem  1  that  Poisson  sampling  is 
alias-free. 

Reconstruction  of  the  correlation  function  C(tr)  of  the  sampled 
process  x(t)  from  the  correlation  sequence  Ci,(n)  is  particularly 


That  Poisson  sampling  is  alias-free  also  follows  from  known  proofs 
of  the  completeness  of  the  Laguerre  functions  without  recourse  to 
Theorem  1. 
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simple  in  this  case,  as  might  be  expected  from  the  fact  that  the 
familiar  Lagiierre  functions  are  obtained  by  orthogonalizing  the 
functions  exp(-f/2),   V  exp{-t/2) ,  Vsxp(-r/2),. . .,  which  differ  only 
trivially  from  (3U).     First  we  recall  the  definition  of  the  Laguerre 
polynomials  (see   [9] j  p.  97) 

(35)  L  (t)   =il^  (tV"^)  =f-  (.l)^(J)^     ,       n>0       . 

n!     dt"  ^  kl 

The  L  (t)  are  orthonormal  over  (0,oo)  with  weight  exp(-t),   i.e. 

^00 

(36)  I     Ljv\{t)exp{-t)dif  =  6^^^     . 

From  (36)  we  obtain  by  a  change  of  variables 

r 

2p  I     L^(2ptr)y2pt)  exp(-2pt»d1^=     6^ 
"P 

so  that  the  functions 

(37)  q^M  =  V^  L^_^(2p*)  exp(-pt)     ,         n  >  1       , 

are  orthonormal  over  (0,oo).  Then,  using  (3U)  and  (35),  we  express 

the  q  (r)  in  terms  of  the  p  (t"),  obtaining 
n  n 

(38)      q^(^)  =  vup  y^  (-i)^(";;^)2^Pk+i^'^^'  ^  -  ^• 

Eq.   (30)   is  the  explicit  form  of  (28)   in  the  case  of  Poisson  sampling. 
Therefore,   in  this  case  (30 )  becomes 

00 
(39)  C(C)     =r^Pnqn^'^^     »      T>0     , 
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where   (see  (29)) 


(w)  p -v^  £ (-1)' h') 2V*i>'    ";i 


Eqs.   (38),   (39)  and  (ItO)   ^ive  0{f)  in  terms  of  the  correlation 
sequence  c,  (n).     The  analogous  expression  for  the  power  spectrum 
F(«)  is  (see   (31)  and  (32)) 

where  the  p     are  given  by  (UO),  and 

n-1 

k  /n-l\  wk+1/ 


T^(a)  =  /57p  Z_^  (-2) 


^]  ^'\o) 


=  >^f^(-2)^f\-^\^p:J''' 


/ .        sn-1 
=  .^ii£l£)_,  n>l         . 

(i^-D)" 

Next  consider  alternate  Poisson  sampling,     where  we  sample  only 
at  every  other  point  generated  by  a  Poisson  process.     In  this  case, 
pit)  is  the  function  Pp(t)  of  (3U),  with  characteristic  function 

0(s)  = P 


(p-is) 

Again,  it  is  easy  to  see  that  this  function  takes  no  value  more  than 
once  on  the  real  axis,  so  that  alternate  Poisson  sampling  is  also 
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alias-free.     liciv'ever,  remarkably  enough,   if  we  s£.r?ple  at  every 
third  time  generated  by  a  Foisson  process,   aliasing  creeps  in.     In 
this  case,   p(1?)  is  the  function  pAv)  of  (3U),  with  characteristic 

function 


(p-is)^ 

Since  this   function  takes  the  same  value  at  the  points  s,   =  p(2  v^  +  i) 
and  s^  =  p(-2/5+i).   Theorem  2  guarantees  the  existence  of  eliases 
is  this  case. 

Let   y  be  the  class  of  probability  densities  in  L     for  which 

pC'V )  =  0    ,  f  <   0    , 

r 

1    tp(t)dr<  00 
o 

(corresponding  to  a  non-zero  rate  of  sampling),   and  whose  Fourier 
^transforms  take  no  value  more  than  once  on  the  real  axis   (and  hence 
are  schlicht  in  the  closed  upper  half-plane).     Each  function  of  (r 
characterizes  an  alias-free  method  of  additive  random  sampling. 
There  sre  vsrious  ways  of  finding  members  of  ^  .     For  example, 
taking  the  real  and  imaginary  parts  of  the  condition 


00  /OO 


'o  'v 


p('C')   exp(i6;)pt)d'C        , 


we  obtain 
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fOO  ,00 


(Ul) 


and 


(U2) 


p(c)cos  co^tdtr    =  p('l?)   cosco^t^'dir 


^00  XX) 


p(ir)  sin  oxjfdf    =     I  p('C)  sin  co^lT  df    . 


If  the  sine  transform  is  positive  for  co  >  C  and  if  the  cosine 
transform  is  strictly  decreasing,  then,  because  of  the  different 
parities  of  the  cosine  and  sine  transforms,    (Ul)  and   0-2)  cannot 
both  hold  for  co^  ^  "2*     Moreover,   if  p(t)  is  strictly  decreasing,   it 
follows  from 

^00  yn/co 


p(f)  sin  co^dr  =         [pCr)  ~p(ir+  -)  +  p(r +—)-...] sin  cotdT 

'o  ^ 


that  its  sine  transform  is  positive  for  oi  >  0.     Thus,   the  class  JP 
of  strictly  decreasing  pCt)  with  strictly  decreasing  cosine  transforms 
is  a  subclass  of  5    ,     Functions  in'^     can  be  found  by  inspecting  tables 
of  Fourier  transforms.     For  example,     p  expC-pt),    (2p/n)exp(-p'^  /n) 
and   (2kp/n)  Bech(kpti)  are  all  in  IP    j  here,  as  usual,   p  is  the  average 
sampling  rate,  and  k  is  the  constant 

"Z^     (2n+l)^ 

Let  p(njt)  ,  n  ■  1,2,...,  be  a  family  of  probability  densities 

in  P  ,  with 
0 

f(jp(njf)df  =   h^ 


-  36  - 

Then  it  is  easily  seen  that  any  function  of  the  form 

00  03 

p('^)   =  V""  a  p(njT)     ,       0  <  a     <  1  ,  V^     a     =  1 

/,         n                            -     n  -  X    .       n 

n=l  n=l 
also  in    a    ,   provided  that 


IS -    - 

o 


,00 

'^  00 


I    Tp(ir)dr  =  ^  a^h^  <   00 
^  n=l 

!fore  generally,   let  pCajtr),  a  >  0,  be  a   family  of  functions  incP  > 
indexed  by  the  continuous  parameter  a,  with 

00 

J    f  p(ajr)df  =  h(a) 

Then  any  function  of  the  form 


j     p(ajf)d5( 


p(r)  =     j     p(ajf  )d$(a)  , 

o 


where  ^(a)  is  a  probability  distribution  function,  is  also  in  q^  , 
provided  that 


,<x> 


'fp('&')d^-  =    h(a)  d}(a)  <  oo 


From  the  foregoing,   it  is   clear  that  alias-free  sampling  methods  exist 
in  abundance. 

We  conclude  by  noting  that  the  rectangular  distribution  with 
probability  density 

(U3)  P^^^  =^     »     Of't:<    2h       , 

p(r)   «0        ,     t'<0,t>2h, 
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and  characteristic  function 

exp(2ihs)  -  1 

(hh)  „  =  j?((s)=  , 

2ihs 

is  not  alias-free.     To  see  this,  we  construct  the  imageJ^in  the 

w-plone  of  the  first  quadrant  of  the  s-plane  under  the  mapping  (hii). 

i±l  is  shown  schematically  in  the  f igurej   it  has  infinitely  many 

disjoint  components  separated  from  one  another  by  the  spiral  curve 

(which  passes  repeatedly  through  the  origin),   and  the  outermost 

component  is  bounded  by  the  interval   (0,1),   the  image  of  the  positive 

imaginary  axis.     Now  consider  two  points  w     and  w_  which,   like  those 

indicated,  lie  in  different  components  of uii.     These  points  correspond 

to  two  points  s,   and  Sp  of  the  half -plane  Ims  >  0,  and  the  straight 

line  segment  joining  s-   and  s_  lies  entirely  in  Ims  >0.     Moreover, 

since  the  image  of  this  segment  must  intersect  the  spiral  at  a  point 

w  ,   there  must  be  a  point  on  the  segment,  s.   say,   at  whJ.ch0(s)  takes 

the  same  value  w-  as  it  does  at  some  real  co     >  0.     Then,   since 
3  o 

w  =  9^(s)  is  analytic  at  s-.  and  therefore  takes  all  values  sufficiently 
close  to  w-  =  0(s-)   in  a  neighborhood  of  s_   (see    [vj ,  p.  18U),   and 
since  0(s)  is  continuous  at  w  ,  we  can  find  some  value  near  w-  which 
is  taken  twice  in  the  open  upper  half-plane.     Hence,  by  Theorem  2, 
the   "  rectangular  sampling  scheme"  with  distribution  (U3)  is  not 
alias-free. 
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0.5 


Image  in  the  w-plane  of  the  first 
qxiadrant  of  the  s-plane  under  the 
mapping   (kh)   with  h  »  1  (schematic). 
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Main  Street 
Haynard,  Mass. 
Attn;  Dr.  John  Ruze 

Radiation,  Inc. 

P.  0.  Drawer  37 

Melbourne ,  Florida 

Attn:  Technical  Library,  Mr.  K.L.  Cox 

Radio  Corp.  of  America 

RCA  Laboratories 

Rocky  Point,  New  York 

Attn:  P.  S.  Carter,  Lab.  Library 

RCA  Laboratories 
David  Sarnoff  Research  Center 
Princeton,  New  Jersey 
Attn:  Miss  ''em  Cloak,  Librarian 
Research  Library 

Radio  Corporation  of  America 
Defense  Electronic  Products 
Building  10,  Floor  7 
Camden  2,  New  Jersey 
Attn:  Mr.  Harold  J.  Schrader 

Staff  Engineer,  Organization 

of  Chief  Technical 

Administrator 

The  Ramo-Wooldridge  Corporation 
P.O.  Box  Ii5!i53  Airport  Station 
Los  Angeles  hS,   California 
Attn:  Margaret  C.  Whitnah, 
Chief  Librarian 

Hoover  Microwave  Co. 
9592  Baltimore  Avenue 
College  Park,  Maryland 

Director,  USAF  Project  RAND 

Via;  Air  Force  Liaison  Office 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 

Rantec  Corporation 

Calabasas,  California 

Attn:  Grace  Keener,  Office  >^anager 

Raytheon  Manufacturing  Company 
Missile  Systems  Division 
Bedford,  Mass. 
Attn:  Mr.  Irving  Goldstein 

Raytheon  Manufacturing  Company 
Wayland  Laboratory,  State  Road 
Wayland,  Mass. 
Attn:  Mr.  Robert  Sorts 

Raytheon  Manufacturing  Comoany 
Wayland  Laboratory 
Wayland,  Mass. 

Attn:  Miss  Alice  0.  Anderson, 
Librarian 

Republic  Aviation  Corporation 
Farmingdale,  Long  Island,  N.  Y. 
Attn:  Engineering  Library 


Thru:  Air  Force  Plant  Representative 
Republic  Aviation  Corp. 
Farmingdale,  Long  Island,  N.Y. 

Rheen  Manufacturing  Company 
9236  East  Hall  Road 
Downey,  California 
Attn:  J.  C.  Joerger 

Trans-Tech,  Inc. 
P.  0.  Box  3h(> 
Frederick,  Maryland 


Ryan  Aeronautical  Company 
Lindbergh  Field 
San  Diego  12,  Calirornia 
Attn:  Librarv  -  unclassified 

Sage  Laboratories 
159  Linden  Street 
Wellesley  81,  Mass. 

Sanders  Associates 
95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  N.  R.  Wild,  Library 

Sandia  Corporation,  Sandia  Base 

P.O.  Box  5900,  Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

Sperry  Gyroscope  Company 

Great  Neck,  Long  Island,  New  York 

Attn:  Florence  W.  Tumbull,  Engr.  Librarian 

Stanford  Research  Institute 

Menlo  Park,  California 

Attn:  Library,  Engineering  Division 

Sylvania  Electric  Products,  Inc. 
100  Plrst  Avenue 
Waltham  5U,  Mass. 

Attn:  Charles  A.  Thornhill,  Report  Librarian 
Waltham  Laboratories  Library 

Systems  Laboratories  Corporation 
liif52  Ventura  Boulevard 
Sherman  Oaks,  California 
Attn:  Donald  L.  Margerum 

TRG,  Inc. 

17  Union  Square  West 

New  York  3,  N.  Y. 

Attn:  M.  L.  Henderson,  Librarian 

A,  S.  Thomas,  Inc. 

I6l  Devonshire  Street 

Boston  10,  Mass, 

Attn;  A.  S.  Thomas,  President 

Bell  Telephone  Laboratories 
Murray  Hill 
New  Jersey 

Chu  Associates 
P.  0.  Box  387 
Whitcomb  Avenue 
Littleton,  Mass. 

Microwave  Associates,  Inc. 
Burlington,  Mass. 

Raytheon  Manufacturing  Company 
Missile  Division 
Hartwell  Road 
Bedford,  Mass. 

Radio  Corporation  of  America 
Aviation  Systems  Laboratory 
225  Crescent  Street 
Waltham,  Mass. 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division  Research  Library 
Box  501j,  Sunnyvale,  California 
Attn:  Miss  Eva  Lou  Robertson, 
Chief  Librarian 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Dr.  W.  C.  Hoffman 

Commander 

AF  Office  of  Scientific  Research 

Air  Research  and  Development  Command 

lUth  Street  and  Constitution  Avenue 

Washington,   D.   C. 

Attn:     Mr.  Otting,   SRY 

Westinghouse  Electric  Corp. 
Electronics  Division 
Friendship  Int'l  Airport  Box  7I46 
Baltimore  3,   Maryland 
Attn:     Engineering  Library 
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Whaeler  Laboratories,  Inc. 
122  Cotter  Hill  Road 
Ireat  Neck,  Nev  York 
Attn:  Mr.  Harold  A.  Wheeler 

Zenith  Plastica  Co. 
Box  91 

Oardena,  California 
Attni  Mr.  S.  S.  Oleesky 

Library  GeophyBical  Institute 
of  the  University  of  Alaska 
College 
Alaska 

Oniversity  of  California 

Berkeley  Ij,  California 

Attn:  Dr.  Samuel  Silver, 

Prof.  Engineering  Science 
Division  of  Elec,  Eng.  Electronics 
Research  Lab. 

nniversity  of  California 
Electronics  Research  Lab, 
332  Cory  Hall 
Berkeley  k,   California 
Attn:  J.  R.  Whiruiery 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
1:800  Oak  Grove  Drive 
Pasadena,  California 
Attn:  Mr.  I.  E.  Mewlan 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
Attn:  Dr.  C.  Papas 

Carnegie  Institute  of  Technology, 
Schenley  Park 

Pittsburgh  13,  Pennsylvania 
Attn:  Prof.  A.  E.  Helns 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,  New  York 

Attn:  Prof.  G.  C.  Dalinan 

University  of  Florida 

Department  of  Electrical  Engineering 

Gainesville,  Florida 

Attn:  Prof.  M.  H.  Latour,  Library 

Library 

Georgia   Institute  of  Technology 

Engineering  Experiment  Station 

Atlanta,    Georgia 

Attn:     Mrs.    J.H.  Crosland,    Librarian 

Harvard  University 

Technical  Reports  Collection 

Gordon  McKay  Library,   303A   Pierce  Hall 

Oxford  Street,   Cambridge  3",   Mass. 

Attn:     Mrs.  E.L.  Hufschraidt,   Librarian 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge  39,   Mass. 
Attn;     Dr.  Fred  L.  Whipple 

University  of   Illinois 
Documents  Division  Library 
Drbana,    Illinois 

University  of  Illinois 
College  of  Engineering 
Urbana,    Illinois 

Attn:     Dr,    P.   E,  Moyes,    Department  of 
Electrical  Engineering 

The  Johns  Hopkins  University 
Horaewood  Campus 
Department  of  Physics 
Baltimore  18,   Maryland 
Attn;     Dr.  Donald  E.  Kerr 

Sandia  Corporation 
Attn:     Organization  l!i23 
Sandia  Base 
AlburjuerTue,   New  Mexico 


Applied  Physics  Laboratory 
The  Johfls Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,   Maryland 
Attn;     Mr.    George  L.   Seielstad 

Massachusetts  Institute  of  Technology 

Research  Laboratory  of  Electronics 

Room  20B-221 

CairtDridge  39,  Massachusetts 

Attn:  John  H.  Hewitt 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

P.  0.  Box  73 

Lexington  73,    Mass. 

Attn:     DociMient  Room  A-229 

University  of  Michigan 
Electronic  Defense  Group 
Engineering  Research  Institute 
Ann  Arbor,   Michigan 
Attn:     J.  A.  Boyd,    Supervisor 

University  of  Michigan 
Engineering  Research  Institute 
Radiation  Laboratory 
Attn:     Prof.   K.   M.   Slegel 
912  N.  Main  St., 
Ann  Arbor,    Michigan 

University  of  Michigan 
Engineering  Research  Institute 
Willow  Run  Laboratories 
Willow  Run  Airport 
Ypsilanti,   Michigan 
Attn:     Librarian 

University  of  Minnesota 

Minneafxjlis  l^,    Minnesota 

Attn:     Mr.   Robert  H.   Stumra,    Library 

Northwestern  University 
Microwave  Laboratories 
Evanston,    Illinois 
Attn:      R.  E.   Bean 

Ohio  State  University  Research  Found. 
Ohio  State  University 
Columbus  10,   Ohio 
Attn:     Dr.  T.E.  Tlce 

Dept.  of  Else.   Engineering 

The  UnlTersity  of  Oklahoma 
Research  Institute 
Norman,   Oklahoma 

Attn;  Prof.  C.  L.  Farrar,  Chairman 
Electrical  Engineering 

Polytechnic  Institute  of  Brooklyn 

Microwave  Research  Institute 

55  Johnson  Street 

Brooklyn,  New  York 

Attn;  Dr.  Arthur  A.  Oliner 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
^^   Johnson  Street 
Brooklyn,  New  York 
Attn;  Hr.  A.  E.  Laommel 

Syracuse  University  Research  Institute 
Collendale  Campus 
Syracuse  10,  New  York 
Attn:  Dr.  C.  S.  Grove,  Jr. 

Director  of  Engineering  Research 

The  University  of  Texas 
Elec.  Engineering  Research  Laboratory 
P.  0.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  Mr.  John  R.  Oerhardt 
Assistant  Director 

The  University  of  Texas 

Defense  Research  Laboratory 

Austin,  Texas 

Attn;  Claude  W.  Horton,  Physics  Library 

Dniverslty  of  Toronto 

Department  of  Electrical  Engineering 

Toronto,  Canada 

Attn;  Prof.  0.  Sinclair 


Lowell  Technological  Institute 
Research  Foundation 
P.  0.  Box  709,  Lowell,  Mass. 
Attn:  Dr.  Charles  R.  Mingins 

University  of  Washington 

Department  of  Electrical  Engineering 

Seattle  5,  Washington 

Attn;  G.  Held,  Associate  Professor 

Stanford  University 
Stanford,  California 
Attn:  Dr.  Chodorow 

Microwave  Laboratory 

Physical  Science  Laboratory 

New  Mexico  College  of  Agriculture 

and  Mechanic  Arts 

State  College,    New  Mexico 

Attn:     Mr.   H.  W.   Haas 

Brown  University 
Department  of  Electrical  Engineering 
Providence,    Rhode  Island 
Attn;     Dr.  C.   M.  Angulo 

Case  Institxjte  of  Technology 

Cleveland,   Ohio 

Attn:      Prof.   S.   Seeley 

Columbia  University 

Department  of  Electrical  Engineering 

Momingside  Heights 

New  York,    N.  Y. 

Attn;  Dr.  Schlesinger 

McGlll  University 

Montreal,  Canada 

Attn;  Prof.  0.  A.  Woonton 

Director,  The  Eaton  Electronics 

Research  Lab. 

Purdue  Dniverslty 

Department  of  Electrical  Engineering 

Lafayette,  Indiana 

Attn:  Dr.  Schulta 

The  Pennsylvania  State  University 
Department  of  Electrical  Engineering 
University  Park,  Pennsylvania 

University  of  Pennsvlvania 

Institute  of  Cooperative  Research 

3!i00  Walnut  Street 

Philadelphia,  Pennsylvania 

Attn:  Dept.  of  Electrical  Engineering 

University  of  Tennessee 
Ferris  Hall 
W.  Cumberland  Avenue 
Knoxvllle  16,  Tennessee 

University  of  Wisconsin 
Department  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn:  Dr.  Scheibe 

University  of  Seattle 

Department  of  Electrical  Engineering 
Seattle,  Washington 
Attn;  Dr.  D.  K.  Reynolds 

Wayne  University 

Detroit,  Michigan 

Attn:  Prof.  A.  F.  Stevenson 

Electronics  Research  Laboratory 
Illinois  Institute  of  Technology 
3300  So.  Federal  Street 
Chicago  16,  Illinois 
Attn:  Dr.  Lester  C.  Peach 
Research  Engineer 

Advisory  Group  on  Electronic  Parts 

Room  103 

Moore  School  Building 

200  South  33rd  Street 

Philadelphia  h,    Pennsylvania 
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Ionosphere  Research  Laboratory 

Pennsylvania  State  College 

State  College,  Pennsylvania 

»TTMi  Professor  A.  H.  Waynick,  Director 


Institute  of  Mathematical  Sciences 

25  Waverly  Place 

New  York  3,  New  York 

ATTNi  Librarian  (3) 

Electronics  Division 
Rand  Corporation 
1700  Main  Street 
Santa  Woniea,  California 
ATTMi  Dr.  Robert  Kalaba 

National  Bureau  of  Standards 

Washington,  D.  C. 

ATTN:  Dr.  W.  K.  Saunders 

Applied  Mathematics  and  Statistics  Lab. 

Stanford  Dnlversity 

Stanford,  (illfomla 

ATTN:  Dr.  Albert  H.  Bowker 

Department  of  Physics  and  Astronomy 
Michigan  State  College 
East  Lansing,  Michigan 
ATTNl  Dr.  A.  Leitner 

Universitv  of  Tennessee 
Knoxvllle,  Tennessee 
ATTN:  Dr.  Frel  A.  Ficken 

Lebanon  Valley  College 
Annvllle,  Pennsylvania 
ATTN:  Professor  B.H.  Bissinger 

General  Atomic 

P.  0.  Box  608 

San  Diego  12     California 

Attn:      Dr.   Edward  Gerjuoy 

Department  of  Physics 
Amherst  College 
Amherst,    Mass. 

ATTN:     Dr.  Arnold  Arons 

California  Institute  of     Technology 
1201  E.  California  Street 
Pasadena,   California 
ATTN:     Dr.  A.  Erdelyl 

Mathematics  Department 

Stanford  University 
Stanford,  California 
ATTN:     Dr.  Harold  Levlne 

University  of  Minnesota 

Minnpapolis  I'j,   Minnesota 

ATTNl     Professor  Paul  C.  Rosenbloom 

Department  of  Mathematics 

Stanford  University 

Stanford,   California 

ATTN:     Professor  Bernard  Epstein 

Applied  Physics  Laboratory 
The  JohnP  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring.   Maryland 
ATTN:     Dr.  B.  S.  Gourary 


(2)Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,   D.  C. 

Electrical  Engineering  Department 
Massachusetts  Institute  of  Technology 
Cambridge  39,   Mass. 
ATTN:     Dr.  L.   J.  Chu 

Nuclear  Developnent  Associates,    Inc. 

5  New  Street 

White  Plains,  New  York 

ATTN:  Library 

California  Institute  of  Technology 
Electrical  Engineering 
Pasadena,   California 
ATTN:     Dr.   Zohrab  A.   Kaprlelian 


Dr.   Rodman  Doll 
311  W.   Cross  Street 
Ypsllanti,   Michigan 

California  Inst,    of  Technology 
Pasadena,   California 
ATTN:     Mr.  Calvin  Wilcox 

Mr.  Robert  Brockhurst 

Woods  Hole  Oceanographio  Institute 

Woods  Hole,   Mass. 

National  Bureau  of  Standards 
Boulder,   Colorado 
ATTN:     Dr.  R.  Gallet 

Dr.  Solomon  L.  Schwebel 

3689  Louis  Road 

Palo  Alto,  California 

University  of  Minnesota 
The  University  of  Library 
Minneapolis  lb,  Minnesota 
ATTN:  Exchange  Division 

Department  of  Mathematics 
University  of  California 
Berkeley,  California 
ATTI:  Professor  Bernard  Friedman 

Lincoln  Laboratory 

Massachusetts  Institute  of  Technology 

P.  0.  Brx  73 

Lexington  73,   Massachusetts 

ATTN:     Dr.  Shou  Chin  Wang,   Room  C-351 

Melpar,  Inc., 

3000  Arlington  Boulevard 

Falls  Church,  Virginia 

ATTN:  Mr.  K.  S.  Kelleher,  Section  Head 

Hq.  Air  Force  Cambridge  Research  Center 

Laure.ice  G.  Hanscom  Field 

Bedford,  Mass. 

.4TTN:  Mr.  Francis  J.  Zucker,  CRRD 


Hq.  Air  Force  Cambridge  Research  Center 

Laurence  G.  Hanscom  Field 

Bedford,  Mass. 

ATTNl   Dr.  Philip  Newman,  CRRK 

Mr.  N.  C.  Gerson 

Trapelo  Road 

South  Lincoln,  Mass. 

Dr.  Richard  B.  Barrar 
Systems  Develo(:mient  Corp. 
2U00  Colorado  Avenue 
Santa  Monica,  Cali_fomia 

Columbia  University  Hudson  Laboratories 
P.O.  Box  239 

lii5  Palisade  Street,   Dobbs  Ferry,   K.  Y. 
ATTNl     Dr.  N.  W.   Johnson 

Institute  of  Fluid  Dynamics 
and  Applied  Mathematics 
University  of  Maryland 
College  Park,   Maryland 
ATTN:     Dr.  Elliott  Montroll 

Department  of  Electrical  Engineering 

Washington  University 

Saint  Louis  5»   Mo. 

ATTN:  Professor  J.  Van  Blade! 

Department  of  the  Navy 

Office  of  Naval  Research  Branch  Office 

1030  E.  Green  Street 

Pasadena  1,  California 

Brandeis  University 
Waltham,  Mass. 
ATTN ;  Library 

General  Electric  Company 
Mictowave  Laboratory 
Electronics  Division 
Stanford  Industrial  Park 
Palo  Alto,  California 
ATTN:  Library 


Smyth  Research  Associates 
3555  Aero  Court 
San  Diego  3,  California 
ATTN:  Dr.  John  B.  Smyth 

Electrical  Engineering 
California  Institute  of  Technology 
Pasadena,  California 
ATTN:  DrJjeorges  G.  Weill 

Naval  Research  Laboratory 

Washington  25,   D.  C. 

ATTN:  Henry  J.  Passerini,  Code  5278A 


Dr.  George  Kear 
5  Culver  Court 
Orinda,   California 

Brooklyn  Polytechnic 

85  Livingston  Street 

Brooklyn,  New  York 

ATTNl  Dr.  Nathan  Marouvltz 

Department  of  Electrical  Engineering 
Brooklyn  Polytechnic 
85  Livingston  Street 
Brooklyn,  New  York 
ATTN:  Dr.  Jerry  Shmoys 

Department  of  Mathematics 
University  of  New  Mexico 
Albuquerque,  New  Mexico 
ATTN:  Dr.  I.  Kolodner 

Mathematics  Department 
Polytechnic  Institute  of  Brooklyn 
Johnson  and  Jay  Street 
Brooklyn,  New  York 
ATTNl  Dr,  Harry  Hochatadt 

Ballistics  Research  Laboratory 
Aberdeen  Proving  Grounds 
Aberdeen,  Maryland 
ATTN:  Dr.  Pullen  Keats 

Dr.  Lester  Kraus 
!j935  Whitehaven  Way 
San  Diego,  California 

University  of  Minnesota 
Institute  of  Technology 
Minneapolis,  Minnesota 
Attn:  Dean  Athelston  Spilhaus 

Ohio  State  University 
Columbus,  Ohio 
Attn:  Prof.  C.  T.  Tad 
Department  of  Electrical  Eng. 

Naval  Research  Laboratories 

Washington  25,  D.  C. 

Attn:  W.  S.  Ament,  Code  5271 

Naval  Research  Laboratory 
Washington  25.  D.  C, 
Attn:  Dr.  Leslie  G.  McCracken,  Jr. 
Code  3933A 

Office  of  Naval  Research 
Department  of  the  Navy 
Attn:  Geophysics  Branch,  Code  ljl6 
Washington  25,  D.  C. 

Office  of  Chief  Signal  Officer 
Signal  Plans  and  Operations  Division 
Attn:  3IGOL-2.  Room  20 
Com.  Liaison  Br.,  Radio  Prop.  Sect. 
Th<?  Pentagon,  Washington  25,  D.  C. 

Defence  Research  Member 
Canadian  Joint  Staff 
2001  Connecticut  Street 
Washington,  D.  C. 

Central  Radio  Prop.  Lab. 
National  Bureau  of  Standards 
Attn:  Technical  Reports  Library 
Boulder,  Colorado 

U.  S.  Weather  Bureau 
U.  S.  Department  of  Commerce 
WashlnfTton  25,  D.  C. 
Attn:  Dr.  Harry  Wexler 
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Federal  Communications  Commission 

Wa«hin,-ton  25,  D.  C. 

Attnt  Mrs.  Barbara  C.  Grimes,  Librarian 

Upper  Atmosphere  Research  Section 
Central  Radio  Propagation  Laboratory 
National  Bureau  of  Standards 
BouMer,  Colorado 

Argonne  National  Laboratory 

P.O.  Box  299 

Lemont,  Illinois 

Attn:  Dr.  Hoylande  D.  Young 


Bell  Telephone  Labs. 
Hurray  Hill,  New  Jersey 
Attn:  Dr.  S.  0.  Rice,  3B 
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Carnegie  Institute  of  Washington 
Dept.  of  Terrestrial  Magnetism 
52U1  Broad  Branch  Road,  N.  W. 
Washington  15,  D.  C. 
Attn:  Library 

Georgia  Tech  Research  Institute 
225  N.  Avenue,  N.  W. 
Attn:  Dr.  James  E.  Boyd 
Atlanta,  Georgia 

University  of  Maryland 
College  Park,  Maryland 
Attnt  Dr.  A.  Weinstein 
Institute  of  Fluid  Dynamics 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

P.  0.  Box  73 

Lexington  73,  Massachusetts 

Attn:  Prof.  Radford,  Division  3  Head 

Willow  Run  Research  Center 
Universitv  of  Michigan 
Willow  Run  Airport 
Tpsilanti,  Michigan 
Attn:  Dr.  C.  L.  Dolph 

School  of  Engineering 
New  York  University 
University  Heights 
New  York,  New  York 

Shell  Fellowship  Committee  of  the 
Shell  Companies  Foundation,  Inc. 
50  West  50th  Street 
New  York  20,  H.  Y. 
Attn:  Mr.  J.  R.  Janssen 

Esso  Research  and  Engineering  Co. 

P.  0.  Box  51 

Linden,  New  Jersey 

Attn:  Mr.  C.  L.  Brown,  Manager 

Union  Carbide  and  Carbon  Corp. 
30  E.  U2nd  Street 
New  York  17,  Now  York 
Attn:  Mr.  L.  E.  Erlandson 

Convair 

San  Diego  12,   California 

Attn:     Mr.  Marvin  Stem 

Bell  Telephone  Labs.,    Inc. 
Ii63  West  Street 
New  York  13,    N.  Y. 
Attn:     Dr.   Hervin  J.   Kelly 

Engineering  Library 
University  of  California 
U05  Hilgard  Avenue 
Los  Angeles  2ii,   California 

Convair,   A  Division  of  General  Dynamics  Corp. 

Daingerfield,   Texas 

Attn:      J.   E.   Arnold,    Division  Manager 

Convair,    A  Division  of  General  Dynamics  Corp. 

San  Diego  12,   California 

Attn:     R.   L.  Bayless,   Chief  Engineer 

Convair,    a  Division  of  General  Dynamics  Corp. 

San  Diego  12,   California 

Attn:     K.    J.  Bossart,   Chief  Engineer-WS107A 


Convair,   A  Div.   of  General  tTnamics 

Corp. 
Fort  Worth  1,   Texas 
Attn:     F.  W.   Davis,   Chief  Engineer 

Convair,.    A  Div.   of  General  Dynamics 

Corp. 
Pomona,   California 
Attn:     C.   D.   Perrine 

Ass 't  Div.   Manager,   Engin. 

Shell  Development  Company 

Exploration  and  Production  Res.   Div. 

3737  Bellalre  Boulevard 

Houston  25,   Texas 

Attn:     Miss  Aphrodite  Mamculldes 

RCA  Laboratories 
Princeton,   New  Jersey 
Attn:     Dr.  Charles  Polk 

Stanford  Research  Institute 
S.   Pasadena,   California 
Attn:     Dr.   J.   Brandstatten 

Wayne  State  University 
Kresge-Hooker  Science  Library 
5250  Second  Boulevard 
Detroit  2,   Michigan 

ARRA 

1  Bond  Street 

Westbury,   L.    I.,    New  York 

Attn:     Dr.   Norman  Spector 

Varlan  Associates 

611  Hansen  Way 

Palo  Alto,   California 

Attn:     Mrs.    Perry  Conway 

Technical  Librarian 

Case   Institute  of  Technology 

Department  of  Electrical  Engin. 

University  Circle 

Cleveland  6,   Ohio 

Attn!     Prof.   Robert   Plonsey 

Dr.  Ming  3.  Wong,    CRRKP 
Air  Force  Cambridge  Research  Center 
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